%% term_grad_U_grad_V
function matrix_A = term_grad_U_grad_V(mesh, fespace)
A1 = term_grad_P_grad_Q(mesh, fespace);
matrix_A = blkdiag(A1, A1);
end
